Comment on "Quantum phase transition in the four-spin exchange antiferromagnet" 
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In a recent paper [Phys. Rev. B80, 174403 (2009)] Kotov et al. studied the paramagnetic- 
to-antiferromagnetic transition in the J-Q model. Their findings were claimed to be in "fairly 
good agreement" with previous quantum Monte-Carlo (QMC) results. In this Comment we show 
that the above claim is misleading and in reality their phase transition point is not only far from 
the corresponding QMC value but also lies in a region of parameter space not yet explored in the 
literature. We also show that their reference dimer state is unstable against formation of a plaquette 
condensate, which could in part explain the large fluctuations they found. 

PACS numbers: 75.10.Jm, 64.70.Tg, 75.40.Cx 



The J-Q model, proposed by Sandviki, is defined by 
the Hamiltonian: 



h = jY,s.s,-qY. Us, - -\ Usi - -\ (1) 




where i, j ... denote sites in a 2D square lattice, N is the 
number of spins, and Si are spin-1/2 operators. The four- 
spin interaction Q and next-nearest neighbor exchange 
J were originally assumed to be positive^. For reasons 
that will become clear later, we will also consider the 
sector with J < (but keeping Q > {)). The model of 
Eq. IT]) with J > exhibits a quantum phase transition 
(QPT) between the antiferromagnetic (AF) and param- 
agnetic singlet phases. However, the location and nature 
of this transition as well as the structure of the singlet 
phase are still debated. In particular, in Ref. Q it was 
concluded, using QMC simulations, that the QPT is 2"^^ 
order, occurs at Qd J ~ 25 and is consistent with the de- 
confined quantum criticality scenario^. In a later work^ 
by the present authors, using a hierarchical mean-field 
(HMF) approach, the location of this QPT was found 
at Qd J ~ 2 — 3, and the transition itself most likely 
describable within the Ginzburg-Landau paradigm. 

The controversy arising from this difference in values 
of Qd J was recently addressed* by Kotov et al.. By con- 
sidering effects of fluctuations around their trial (param- 
agnetic) columnar dimer state (CDS), they came to the 
conclusion that a mean-field theory is incapable of cor- 
rectly describing the QPT in the J-Q model. The nature 
and numerical value of the phase transition point, found 
in their paper, was claimed to agree with the work of 
Sandviki. This circumstance was also used to speculate 
that the small, compared to QMC, value of Qd J result- 
ing from the analysis of Ref. 0, is due to limitations of 
the HMF approach. We would like to stress that the 
HMF method is not a standard mean-field approach, as 
the one used in Ref. 4, and it becomes an exact method 
in the thermodynamic limit. Moreover, for a finite sys- 
tem HMF can be implemented as a variational theory in 
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Figure 1. Qc/J (Qc > 0), obtained in Refs. audi. 

terms of the energy, and finite-size scaling needs to be 
performed to extrapolate to the thermodynamic limit. 

In the present note we show that the claims of Ref. 
are misleading. We strongly oppose the statement 
"Near the QCP [quantum critical point], whose location 
[Kd Jk ~ 2.16] we find in fairly good agreement with re- 
cent QMC studies...", made in Ref. 0, by demonstrating 
that a direct comparison of the two results is inappropri- 
ate, because the phase transition point, claimed by Kotov 
et al., in reality lies in the yet unexplored region J < 
of the J-Q phase diagram. This circumstance may raise 
doubts regarding the nature of the magnetic state found 
by Kotov et al. We use a simple variational argument 
to show that the large fluctuations found in their work 
could be attributed to an instability of their reference 
dimer state against formation of a plaquette condensate. 

Up to an irrelevant constant, the Hamiltonian ([1]) can 
be rewritten* in the form: 

Hk = JkY.^^,)S^S, ~ KY.(,,u) {S^S,){SuSi) (2) 

with Jk, K > Q. The new coupling constant K/Jk is 
related to the old one Q/J by the formula: 

Q/J - {K/Jk)/[1 - K/2Jk\ (3) 

Clearly, for positive Q and J the parameter range that 
can be explored with the Hamiltonian is ^ K/ Jk ^ 
2. The values of Kj Jk, larger than 2 correspond to the 
region Qj J < 0, which implies either (i) Q < 0, J > 0, 
or (ii) Q > 0, J < 0. We will assume, as Ref. 11, that K 
and Q have the same sign and disregard the case (i). The 
case (ii) defines the ferromagnetic (FM) part of the phase 
diagram of the Hamiltonian ([T]). In the new representa- 
tion ([3]) the results of Refs. [iTandl^ are Kd Jk ~ 1-85 
and Kc/Jk ~ 1 — 1.2, respectively. The critical value. 
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Figure 2. Phase diagram of the J-Q model for Q > 0. The 
part J > was computed in Ref. [3- The point Qo/ J « —1.3 
indicates a T'' order QPT between the FM state and singlet 
plaquette crystal (PC) phase'^. The arrow shows Qc/ J ~ —27, 
obtained in Ref. jj. The PC regions at J < and J > are 
adiabatically connected by changing J {Q/J) through zero 
(infinity), i.e. there is no direct FM-to-AF transition. 



obtained in Ref. is Kc/ Jk ~ 2.16 which, after going 
back to the original units, corresponds to Qd J ~ —27. 
This value should be compared with the result of QMC 
simulations^ Qd J ~ +25. Results of Refs. [H, [3 and|3 
are summarized in Fig. [T] Thus, the claimi that the 
transition point Kc/ Jk ~ 2.16 is in "good agreement" 
with the QMC result, is unjustified. A fair comparison 
involves more than numerology, as we show below. 



Table I. Numerical values for the l'^' order transition point 
Qo/J, obtained using ED and HMF for different cluster sizes. 
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Figure 3. Energy ([SJ as a function of d. The plaquette state 
corresponds to a local minimum £o(57r/6) = —73/96. The 
columnar dimer configuration has a higher energy £o(0) — 
—66/96 and does not correspond to an extremal point of eo{d). 



Since the ratio Qc/ J obtained by Kotov et al. is neg- 
ative, one might think that in reality Ref. 4 studies the 
FM-to-singlet phase QPT. However, this is not necessar- 
ily the case. In fact, negative values of Q/J [Q > 0) do 
not imply a FM phase. In order to address this issue, 
we used the HMF approach'^ and exact diagonalization 
(ED), both in clusters of 2 x 2 and 4x4 spins, to de- 
termine the J-Q phase diagram in the AF [J > 0) and 
FM (J < 0) regimes. The results are shown in Fig. [5] 
The QPT separating FM and singlet phases is I*'* or- 
der: it manifests itself as a level crossing both in ED and 
HMF. In Tabic H] we present numerical values for Qo/J 
obtained from these methods. We see that they are in 
excellent agreement with each other, which is not sur- 
prising, given the fact that the FM state is semiclassical 
and the singlet phase is gapped. Thus, the paramagnetic 
phase displays two, FM and AF instabilities. The region 
K/Jk > 2 and the point Kc/Jk ~ 2.16 reside in the 
singlet phase of Fig. [21 at negative values of Q/J. On 
the other hand, the CDS^ can accomodate a Neel phase 
only if the magnetic unit cell includes two dimers. Ref. 
3 seems to consider only homogeneous phases with one 
dimer per unit cell. This might raise doubts regarding 
the AF nature of the CDS instability discussed by Kotov 
et al. 



Finally, we shall comment on the huge fluctuation cor- 



rections to the value of Kc, found in Ref. |3- These 
corrections were argued to be responsible for shifting the 
QPT towards negative values of Q/J and seem to be in- 
timately related to the use by Kotov et al. of the CDS 
as a physical vacuum for their analysis. This dimer state 
is unstable, when compared to the plaquette structure^. 
Let us consider a trial paramagnetic state of the form: 

l*o) = l[[cos e\^M)+ sin d{\ija) - {■'PdHa)^], (4) 
□ 

where the product runs over N/A 2x2 plaquettes. The 
wavefunction of Eq. ^ interpolates between the dimer 
condensate^- Yl \ipd) (iV'd) is the direct product of spin 
singlets on two parallel links of a plaquette), for 6* = 0, 
and the plaquette stated DlV'n); which corresponds to 
6 = 57r/6. The expectation value of Hamiltonian ([Ij, 
eo{e) = {^omm-^om/NQ, is given by: 

ea{0) = -{J/8Q) {sine - VS cos ey - (l/96)x (5) 

X (51 + 13 cos 26* + 2 cos 40- 13\/3sin26' + 2\/3 sin 49). 

For Q ^ J, this function is shown in Fig. [3] The plaque- 
tte state corresponds to a local energy minimum. On the 
contrary, the CDS does not describe any extremal point. 
Moreover, it has a higher energy compared to the pla- 
quette configuration. The energy difference between the 
CDS and a correlated plaquette structure will be even 
larger if 4 x 4 spin clusters^, are used as a basis for the 
HMF analysis. 
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